Energy Levels of Interacting Fields in a Box 
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We study the influence of boundary conditions on energy levels of interacting fields 
in a box and discuss some consequences when we change the size of the box. In 
order to do this we calculate the energy levels of bound states of a scalar massive 
field x interacting with another scalar field <f> through the lagrangian C int = \g4> 2 X 2 
in an one-dimensional box, on which we impose Dirichlet boundary conditions. We 
have found that the gap between the bound states changes with the size of the box 
in a non-trivial way. For the case the masses of the two fields are equal and for 
large box the energy levels of Dashen-Hasslacher-Neveu (DHN model) (Dashen et 
al, 1974) are recovered and we have a kind of boson condensate for the ground state. 
Below to a critical box size L ~ 2.93 the ground state level splits, which we 
interpret as particle-antiparticle production under small perturbations of box size. 
Below another critical sizes (L ~ an d (L ~ 1.71^?) of the box, the ground 

state and first excited state merge in the continuum part of the spectrum. 
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1. INTRODUCTION 



Consider a simple system of two interacting fields, described by the lagrangian density 



where A and g are coupling contants. In this work we study the influence of boundary conditions 
on energy levels of the field x taking into account its interaction with field <fi. From above la- 
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grangian we derive the following equations of motion for both fields. 

-a^ X + M^-3# 2 x = 0, (1) 

-a^0 + M|0-A0 3 = O. (2) 

In Eq.(2) we have neglected the term 3 g<px 2 which can be interpreted as the back-reaction 
of field x on the mass term of <fi. This can be achieved if we impose that |%| << min{-^^}. 



Of course other regimes can be studied from Eq. (§) adopting different approximations. 

It is well known the static solutions of the classical equation of motion (fj) are given by 
(Dashen et al, 1974) 

These are Kink solutions which connect two vacua at x — ±oo. In this work unlike x, the 
field (f> is not changed by boundary conditions. So in our approximation they are transparent 
to (p. In a future work we will consider the most general case. 

Substituting these solutions in Eq.(P we obtain 

-d»d, X + Mix ~ ^M 2 tanh 2 (^-) X = 0. 

Since we are interested in stationary solutions, we can write xi x ) = e~ tuJt/ if)(x), where u) are 
energy eigenvalues and the previous equation reduces to 

^(x) + (m 2 + u, 2 - 3 |M, 2 tanh\^L)^ tP(x) = 0. (3) 

This equation is similar to that one from DHN model (Dashen et al, 1974), which describe 
Kinks in (1 + 1) dimensions, but here two different mass parameters appear in the potential. 

In the section 2, we calculate the bounded energy levels of the field x ( or VO constrained 
to a box in (1 + 1) dimension with Dirichlet boundary conditions. We discuss how we can 
interpret the splitting of the ground state when the box is shortened below a critical value, as 
particle production from a vacuum condensate. 

2. BOUNDED SPECTRUM 



In this section we obtain the energy levels of field x imposing Dirichlet boundary conditions. 

V2 



This is as follows. With the changing of variables z = ^4=- and uo 2 = — ^— - the equation 



can be written as 



ip(z) + + e - 2 - Qhanh 2 z\ ip(z) 
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Making a new changing of variable, namely, E = 2(3 + e — 2, where we have defined the 

M x 



M 2 

mass ratio (3 = -rj£, the above equation reduces to 



— V>0) + (£ - 6^-tanh 2 zU(z) = 0. (4) 
az 2 A 

For discrete (bounded) spectrum case, < uj 2 < 2M 2 , we have that 2(3 < E < 6(3. 
In order to find the corresponding solution of the equation (U) we make the following variable 
dependent transformation (Morse and Fesbach, 1953): 

ip(z) = sech k (z)Y(z), (5) 

where the parameter k G R will be determined below. 
Substituting (||) in (£|) we get an equation for Y(z) 

d 2 Y „ , , , , dY 



— 2k tanh(z)— — h 



dz 2 dz 



k 2 + E-6{ | 6 g fc2 fc 
sech 2 (z) X 



sech 2 {z)Y(z) = 0. (6) 



In this work, we discuss the particular case g = A. The general case, including an asymptotic 
study for weak and strong coupling constant g will be presented elsewhere. Nevertheless the 
above condition leads to interesting results. Firstly it is possible to turn this equation into a 
hypergeometric differential equation. This is as follows: 

1. - We impose the term dependent on the variable z in square brackets, equal to zero, which 
gives a relation between k and E, i.e., 

k = ±V6 - E. (7) 

2. - Making a changing of the independent variable, namely, 

1 

(i= -(1 -tanh(z)) (8) 

we obtain a hypergeometrical differential equation 

d 2 Y dY 
H(l - n)—^ + (k + 1 — 2(k + l)//)-r- - (k + 3)(k - 2)Y(fi) = 0. (9) 

This equation has regular singular points at ji = 0, fi = 1 and /j, — oo, and parameters 
a = & + 3, b = k — 2 and c = fc + 1. The two independent analytic solutions are given by 
(Gradshteyn and Ryzhik, 1980): 

Y 1 { f i)= 2 F 1 {k + 3,k-2;k + l]n), (10) 
Y 2 (fi) =ti~ k 2 F 1 (-2,3; 1- k; fj). (11) 
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As it is well known, if the parameter c is a positive integer a solution of the equation (^|) will 
be Yi(n) and, if c is a negative integer, it will be given by Y 2 (n). Also there are other solutions 
containing a logarithmic term (Fuchs Theorem). If c is a non-integer number, the set of two 
solutions Yi(fi) and y 2 (/i) is a system of linearly independent solutions (Butkov, 1968). Next 
we study these three different cases for c. From now on, in order to simplify the notation, we 
denote the hypergeometric function 2 Fi simply as F. 

CASE I: c Is a Positive Integer. 

In this case, c = k + 1 = n, where n = 1, 2, 3, . . ., and by Fuchs Theorem (Butkov, 1968), 
the general solution is given by 

oo 

Y{/j) = AYy(n) + BY x (n) In | // | +£ £ a s (-k)/j 1 s ' k , 

s=0 

where A and B are arbitrary constants. 

It is easy to see that this solution has a logarithmic divergence at [i — 0. Since we want 
the solution Y{p) be bounded, we must take B = 0. This comes from the constraint that for a 
large box (L — ► oo) we must recover DHN's Kink and its bounded energy levels. So the solution 
reduces to 

Y((j) = AF(n + 2,n-3;n;fi), 

where k was substituted by n — 1. 

Using the relations @ and (0), as well changing of variable z = in the above expression, 
we obtain 

V>(x) = A sech {n - 1 \^-) F{n + 2, n - 3; n; i{l - tanh(^)}). (12) 

Below, in order to simplify the notation, we will denote M<f, as M. 

Now we impose Dirichlet boundary conditions at x = ±-| for the solution ([12]), that is 

tf(T§) = Aaech^C^) F(n + 2, n - 3; n; i{l ± tonft(^)}) = 0. 

From these relations we get the condition 

1- , , ML . , , . , 

F(n + 2, n - 3; n; -{1 ± tan/i(^=)}) = 0. (13) 

On the other hand, since 2/3 < E < 6/5, from equation (H), we can show that the parameter 
fc satisfies the inequalities 



6(1 - p) < k < J6 - 2(3 or - JG - 2(3 < k < -v/6(l - /?). 
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Since we have considered that k G R, then from the above relations we obtain f3 < 1. 

So, the possible values of parameter k are in the intervals < k < v6 ~ 2, 44 or —2, 44 ~ 
— v6 < A; < 0. As in this case c is a positive integer, we must take the interval < k < \^6. 

The allowed integer values of k in the interval < k < \/6 and the correspondent values of 
E (from Eq.([I3J)) are (using the relation k — n — 1) 



for n = 2, k = 1 then E = 5 



for n = 3, k = 2 then E = 2. (14) 



The hypergeometric function in ( |13l) can be written as Jacobi Polinomials (Abramowitz and 
Stegun, 1972) and it is not difficult to prove that substituting the allowed values of k ( |I~4"D into 
(P~3|) we do not get any consistent solution. 

CASE II: c Is a Negative Integer. 

Now we consider c = k + 1 = — n, where n = 1,2, 3,..., and then we must consider the 
interval — \/6 < k < 0. The solution of the hypergeometric equation, in this case, is given by 
Y2 (yu). As for the previous case, by Fuchs Theorem, we have that the general solution is given 
by 



Y(fj) = AY 2 (fJ.) +BY 2 (n)]n | /1 | +B £ a s (-k)fi s -\ 

s=0 

where A and B are arbitrary constants. 

In this case we can not discard the S-terms in the same way as in the previous case, since 
for /i — > we have n~ k In /i — > and s — k > for s = 0, 1, 2, • • • . So we don't have a singularity 
at fj, — 0. Nevertheless, we notice that the relation is divergent in the asymptotic limit of 
fi = 0, i.e. 

hm . , -» 00. 

m-o y 2 (^) 

Again, we can impose our natural boundary condition, that is, for very large box (L — > 00), 
the DHN solution (Dashen et al, 1974) should be recovered. In order to do this, we must impose 
the following asymptotic condition 

H7 2 (/i) 

In order to this condition be valid, the coefficient B must vanish. Therefore, our solution 
reduces to 

F(/i) = A^ n+ V F(-2, 3; 2 + n; fi), 
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where k was substituted by — (n + 1), and A ^ 0. 

As we did before, using the relations (§) and (|^), as well the changing of variable z = y| 
in the above expression, we can write the solution in the original variables: 

, s Mx (1 -tanh(^£)Y n+ V i Mx 

^{x) = Asech~^\^) X - F(-2, 3; 2 + n; -{1 - tan/>(^)}). (15) 

Now, imposing Dirichlet boundary conditions at x = ±|f, namely 

L , ML {1 ±tan/i(^)}(" +1 ) i ML 

^ v 2 ; v 2v^ S^ 1 ) v ' ' '2 l y 2^2 

we obtain the condition 

F(-2, 3; 2 + n; ^{1 ± tanh(^)}) = 0. (16) 

In this case, the allowed integer values of k, in the interval — \/6 < k < are using the 
relation k = —(n + 1) given by 

for n = 0, k = — 1 then E = 5, 



for n = 1, k = -2 then E = 2. (17) 



As in the previous case it is not difficult to write the above hypergeometric function ( ]T6| ) as 
Jacobi Polinomials (Abramowitz and Stegun, 1972). A quick analysis shows that no consistent 
solution exists for finite L ^ 0. For L = oo, the DHN's case (n = 0) is obtained. 

CASE III: c Is a Non Integer. 

Since c is a non- integer number (positive or negative), the general solution is given by 

F(/i) = A F(k + 3, k - 2; k + l; fi) + B (fi)- k F(-2,3;l - k; fi). (18) 

Now we repeat the previous steps, i.e., using the relations @ and (^), as well the change 
of variable z = ^= in the previous expression, and obtain 

ip(x) = sech k (^) (a F(k + 3, k - 2; k + 1; ~{1 - tanh(^)}) + 



Mr 1 Mr \ 

B 2 k [l - tanh(-^)}- k F(-2, 3; 1 - *; -{1 - tanh(-^)})\ . (19) 
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Now must determine k and, afterwards, using relation (|7|), find the allowed values of E. 
The Dirichlet boundary conditions at x — ±~ are given by 

^(=f| ) = M- F(k + 3, fc - 2; fc + 1; ^{1 ± tanh(^=)}) + 

5 2 k [l ± tanM^|)]- fe F(-2, 3; 1 - *; ±{1 ± tanh(^L)})\ = 0. 
From these relations we obtain 

AF(k + 3,k — 2;k + l; ^{1 ± *arafc(^|)}) + 

5 2 fc [l ± tanM^|)]- fc F(-2, 3; 1 - *; ±{1 ± tanh(^L)}) = 0. 

This is a system of homogeneous equations for A and B. So this system has a non-trivial 
solution only if the determinant of system is zero, i.e., 

{1 - tanh(^)}- h F(k + 3,k-2;k + 1; ~{1 + tanh(—^)}) x 
F(-2, 3; 1 - *; i{l - tanh(^=)}) - {1 + tanA(^)}-*x 

F(fc + 3, k - 2; + 1; ~{1 - tanh(^L)}) F(-2, 3; 1 - fc; ~{1 + tan/i(^|)}) = 0. (20) 

Again we can write the above hypergeometric functions as Jacobi Polinomials (Abramowitz 
and Stegun, 1972). In this way, from (|20|) we obtain a transcendental equation for the parameter 
k, i.e., 

+ tanh(M±) \ k ( k*-l + 3ktanh(ML) + 3tanh*(M%) \ 
l-tanh($%)) \k 2 -l-3ktanh(M±)+3tanh 2 (M±)) ' 1 J 

Notice that substituting k by —A; in equation (pT]), it can be verified that this equation 
is satisfied. So these solutions are valid in the intervals (0,1) U (1,2) U (2,^/6) C R and 
(— \^6, —2) U (—2, —1) U (—1, 0) C R. Observe that for the negative sign, k = — 1 is a solution 
of equation (|2l|) , but it is not allowed since now we are considering k Z. 

An analytic solution for the transcendental equation (|2l|) was not founded but it is possible 
to obtain numerical solutions for it as we can see in the figure below for the case that (3 = 1. 
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Fig. 1. Shifts of Bound States with the size of the Box. Here (3 = 1. The 
horizontal lines show the asymptotic values of the DHN model (Dashen 
et al, 1974). 



The Fig. 1 shows the relation between parameter k and the size L of the box. Observe 
that both ground state and the first excited state shift with the size of the box. Starting 
with a very large box size L, as it decreases, the values of c^o and u>\ increases from Uq = 
(k = —2 and L — > oo) to u>q = \[2M (k = 0), and from u% = \J~\M (k = —1, and L — > oo) 
to uj\ = \plM (k = 0). Close to critical value ~ 0.6 for u)q, and ~ 1.71 for u)\ these bound 
states merge in the continuum part of the spectrum (k = 0). On the other hand, for L — > oo 
the values of ujq and u\ decreases until to reach their minimum values ujq = and u)\ = J\M 
respectively, just as in the DHN model. This behaviour happens for both cases of positive and 
negative signs in the equation (^1|) with negative k. For any positive k and large box sizes L 



the equation (£]]) has no solution, which can be checked analitically (see also the Fig.l above). 
So we discard positive values of k since for a large box (L — > oo) they do not lead to DHN's 
energy levels. 

We can also obtain a relation between ^jp- and box size L. From equation (0) and the 
relation oj 1 = ^ £ ~ 2 ^ M we get 
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In the Fig. 2 below we plotted two cases for this relation, namely (3—1 and (3 — \ . 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 
ML/2.82 

Fig. 2. Energy levels of the ground state of field x ( or VO dependent of 
the box size L. The horizontal lines show the asymptotic values of DHN 
model (Dashen et al, 1974). 

Now we can see that for the first case ((3 = 1), the ground state u)q, for large box size may 
form an aproximated vacuum particle-antiparticle condensate, since both levels are exponen- 
tially close one from another. Roughly speaking below a critical size (~ 2.93) any perturbation 
of box size will induce pair formation from this energy level. This could be important for 
particle production in the presence of strong fields. Squeezed fields could have an enhanced 
production. 

Another interesting behaviour of the energy levels is showed in Fig. 3 below. We plotted 
the difference between the levels u>q and oui, for arbitrary distance L. The interesting fact about 
this is the peak around the critical value ~ 2.93 and that the increasing part of the curve 
(left-hand side) from the peak shows a not smooth, growth with several secondary maxima and 
minima. 
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Fig. 3. The gap between levels u and u\ versus the size L of the box. 

3. CONCLUSIONS 

In this work we calculated the solutions of a Klein-Gordon type equation in an one- dimensional 
box for which we impose Dirichlet boundary conditions, and in the presence of a Kink-type 
potential generated by a second scalar, self-interacting field which, in our approximation, is not 
subjected to the boundary conditions. Energy levels of bounded states for non-trivial solutions 
are obtained as roots of a transcendental equation involving L and k. Although we have not 
obtained analytic solutions to it, we have studied numerical solutions (see Fig. 1). 

The ground state ujq and the first excited state uj\ shifts with the size L of the box (see fig. 
1), for both cases of positive and negative signs. As the size L of the box decreases, uq increases 
in the interval [0, V2M], and uj\ in [^/|M, \p2M\. Close to the critical value ~ 0.6 for uq, and 
~ 1.71 for u\, all the bounded states merge in the continuum part of the spectrum. For large 
distances (L = oo) we obtain the energy levels of DHN model. The decrease of L induces shifts 
on the bound states levels of the system, and close to a critical size ~ 2.93 we have "just barely 
bound" condensate (Morse and Feshbach, 1953) that may decay against a small perturbation 
on the system, with particle pair creation (Fig. 2). 

The gap between the two bound states presents a peak at ~ 2.93 and shows a wildly 
(non-smooth) behaviour. It is interesting that the critical value for the splitting of the levels 



Aco(L)=co 1 (L)-co (L) 
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coincide with the value to the peak position, but we do not have any explanation for this fact 
up to now. This, as well the study of the system taking into account finite boundary conditions 
to both fields (ft and x will be done elsewhere. 
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